Screened hydro-dynamic interaction in a narrow channel 
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We study experimentally and theoretically the hydro-dynamic coupling between Brownian colloidal 
particles diffusing along a linear channel. The quasi-one-dimensional confinement, unlike other 
constrained geometries, leads to a sharply screened interaction. Consequently, particles move in 
concert only when their mutual distance is smaller than the channel width, and two-body interactions 
remain dominant up to high particle densities. The coupling in a cylindrical channel is predicted to 
reverse sign at a certain distance, yet this unusual effect is too small to be currently detectable. 

PACS numbers: 83.50.Ha, 83.80.Hj, 82.70.Dd 



Diffusion along narrow channels is encountered in vari- 
ous circumstances, such as transport in porous materials 
jl| and penetration through biological ion channels ||. 
Much attention has been devoted to the correlated mo- 
tion in so-called single-file systems due to the inability 
of particles to bypass one another, leading to anomalous 
diffusion jl], [jj. The anomalous regime, however, sets in 
only at long enough times when particle collisions become 
appreciable Q. Yet, if the channel is filled with liquid, 
the motion of otherwise non-interacting Brownian par- 
ticles can become correlated at times shorter than the 
collision time through the flow field that they create. In- 
deed, when we observed the diffusion of micron-size col- 
loidal particles in a water-filled channel (Fig. [l]), the most 
striking feature was the occasional concerted motion of 
several close-by particles in a dynamical "train" lasting 
for a few seconds. In the current Letter we present a 
quantitative study of this coupling. 

When a particle moves through a fluid it creates a flow 
that affects the velocities of other particles in its vicin- 
ity. Recently there has been renewed interest in such 
hydrodynamic interactions, particularly in confined ge- 
ometries, due to their role in the behavior of colloidal 
suspensions and the development of techniques using dig- 
ital video microscopy 

H I; ! ! 0- Colloidal parti- 
cles in a finite container H, near a single wall ||, @, D , 
and between two walls Jq, |9|] were studied. The hydro- 
dynamic interactions in those geometries are always at- 
tractive (i.e., creating positive velocity correlations) and 
long-ranged: in an unbounded fluid the interaction de- 
cays with inter-particle distance x as 1/x pi) ; the one 
between particles moving near and parallel to a single 
wall decays as 1/x 3 |?], |l2); and for particles moving be- 
tween and parallel to two walls it decays as 1/x 2 [ p"3|| . 
More constrained geometries — perpendicular to the walls 
in a two-wall configuration |l3| and along a cylindrical 
tube jl4], [l5| — are essentially different, in that point dis- 
turbances should create flows with an exponential spatial 
decay. Here we demonstrate that confinement in a linear 
channel indeed sharply screens the hydrodynamic inter- 
action and may even change its sign. 



The experimental system consists of an aqueous sus- 
pension of silica colloidal spheres (density 2.2 g/cm 3 ) con- 
fined in long narrow grooves (see Fig. [j] and inset of Fig. 
|]). The grooves were printed on a polydimethysiloxane 
substrate from a master pattern fabricated lithographi- 
cally on a Si wafer (Stanford Nanofabrication Facility). 
A drop of suspension was enclosed between the polymer 
mold and a cover slip with a spacer (~100 /xm), so that 
the top of the groove was open to a layer of fluid. Digi- 
tal video microscopy was used to extract time-dependent 
two-dimensional trajectories of the spheres (time resolu- 
tion 0.033 s). Details of sample preparation and data 
analysis were described elsewhere [0, Eg. 
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FIG. 1: Microscope images of silica colloidal spheres confined 
in linear grooves (cf. Table [|): a) sample 1, T) = 0.17; b) 
sample 2, n — 0.35; c) sample 3, n — 0.36. 



no. 


a ((im) 


w (/jm) 


h (fim) 


L (mm) 


1 


0.79 ±0.02 


3.0 ±0.3 


3.0 ±0.3 


2 


2 


1.85 ±0.05 


5.0 ±0.1 


4.0 ±0.5 


10 


3 


2.57 ±0.1 


7.0 ±0.1 


4.0 ±0.5 


10 



TABLE I: Parameters of the three experimental samples. 

We studied samples of three different parameter sets, 
i.e., different values of sphere radius a, channel width w, 
channel depth h, and channel length L; see Table | and 
Fig. [I]. (Spheres of sample 1 were manufactured by Duke 
Scientific and those of samples 2 and 3 by Bangs Labora- 
tory.) The ratio a/w ranges between 0.26 and 0.37; it is 
always larger than 1/4, such that spheres cannot bypass 
one another. (Higher ratios pose a difficulty in loading 
the spheres into the channel.) The linear packing fraction 
r\ was determined as r\ = 2Na/l, where I is the length of 
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the channel section in the field of view (I = 106/im for 
sample 1 and 220/xm for samples 2 and 3), and N is the 
number of spheres in that section. 

Studies of the equilibrium structure and single-particle 
dynamics of sample 1 were previously reported (J, [l6). 
These experiments show that the colloidal motion is 
tightly confined to the center of the groove, with trans- 
verse fluctuations of less than 0.2a. The particle pair- 
potential consists of a short-ranged (screened electro- 
static) repulsion, extending to surface-surface separa- 
tions of ~ 0.2 /xm, followed by a weak attractive well 
of ~ — O.SfceT 1 (ksT being the thermal energy) at a sepa- 
ration of ~ 0.3 yitm. At the range of the dynamic coupling 
reported here the pair-potential is practically zero |L6| . 
We thus neglect in the current study any transverse mo- 
tion and direct interactions. 

We choose to characterize the coupling by the collective 
and relative diffusion of particle pairs, i.e., the fluctua- 
tions of their center of mass and mutual distance. We 
define the collective diffusivity, D + , and the relative one, 
D~ , as 

D ± (x) = ([Ax 1 (t)±Ax 2 (t)] 2 )/U, (1) 

where Axi is the displacement of particle i during time t 
and x is the short-time average distance between the cen- 
ters of the two spheres. The coefficients have been defined 
such that, in the absence of coupling, they both reduce 
to the self-diffusivity of a single particle, D s . In our data 
analysis D ± (x) were calculated from the histogram of 
short-time (t < 0.2 s) trajectories, i.e., the slope of the 
mean-square displacements, [l/4N(x)] EiL^ [ Aa; i (*) ± 
Axi-i(t)] 2 , as a function of t, where N(x) is the number 
of nearest-neighbor pairs whose mutual distance x falls 
in the range (x — 8x/2, x + Sx/2). (We took Sx = 0.22a 
for all samples.) The tracking time t must be kept suf- 
ficiently short, so that the spheres do not cross over to 
the anomalous subdiffusion regime ji) , and their diffusion 
distance remains smaller than Sx, so x can be assumed 
constant. We verified that, within experimental error, 
the measured slopes were constant at least in the range 
t = 0.1-0.5 s. 

Figure ^ shows the measured collective and relative 
diffusivities for the various samples. The diffusivities 
sharply decay to D s for inter-particle distances larger 
than about twice the channel width. Moreover, the cou- 
pling is practically insensitive to changes in density up 
to high values of 77, implying that three-body and higher 
terms remain negligible. At r\ = 0.61 sample 1 exhibits a 
considerable shift of the entire curves towards lower diffu- 
sivity, i.e., a decrease in D s |t7|. These observations are 
in essential contrast with the case of less confined geome- 
tries, where the hydrodynamic interaction is long-ranged 
with significant many-body effects 

We now turn to a theoretical estimate of the hydro- 
dynamic coupling. The Reynolds number of the system 
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FIG. 2: Collective (D + , upper curves) and relative (D~ , lower 
curves) pair diffusivities as a function of inter-particle distance 
x for various values of packing fraction r\. Figures a, b and c 
correspond to samples 1, 2 and 3, respectively. 

is of order 10~ 6 . We need, in principle, to calculate the 
Stokes flow due to the motion of two particles, subject to 
no-slip boundary conditions at the surfaces of the chan- 
nel and particles Jul |l8). This is technically very hard 
and, furthermore, we seek a description of more general 
applicability. We resort to three simplifications: (i) The 
actual geometry of a partially open, rectangular channel 
is replaced by an effective cylinder, of diameter 2R = 710, 
where 7 is a geometrical prefactor of order 1 , to be treated 
as a fitting parameter. We expect to obtain one value of 
7 for a square cross-section (samples 1 and 2), and an- 
other for a rectangular one (sample 3). Because of the 
open top of the actual channel, we also expect the effec- 
tive cylinder to be significantly wider, i.e., 7 > 1. (ii) 
The particle motion is assumed to be restricted to the 
central axis of the channel. (As has been noted above, 
this is a good approximation in our case.) (iii) The par- 
ticle size is assumed to be much smaller than both the 
channel width, a <C R, and the inter-particle distance, 
a <C x. This requirement, though not strictly fulfilled in 
practice, allows us to treat the effect of one particle on 
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the flow near the other and near the walls as if it were 
exerting a point force on the fluid, i.e., to consider merely 
the fundamental solution (Green's function) of the Stokes 
flow in the channel. (This is sometimes referred to as the 
stokeslet approximation p"2|.) 

The displacement fluctuations of a particle pair 
are used to define a two-particle diffusion tensor, 
(Axi(t)Axj(t)) = 2Dij(x)t, i,j — 1,2, which can be de- 
composed into a self-diffusion term and a coupling one: 
Dij = SijD s + (1 — Sij)D c (x). The two eigenvalues, 
D ± = D s ±D c (x), are the collective and relative diffusiv- 
ities of the pair as defined in Eq. (Q) . Usually particles 
entrain one another and D c > 0, i.e., the hydrodynamic 
interaction enhances the collective mode and suppresses 
the relative one. 

To leading order in a/ R and a/x, the flow field induced 
by the motion of one particle in the vicinity of the other 
is that of a point force. This fluid velocity per unit force 
is the change in the mobility of the entrained particle. 
Thus D c (x) ~ k^TG{x), where G is the xx component 
of the Stokes-flow Oseen tensor (Green's function) in a 
channel. Using the diffusivity in an unbounded fluid, 
Do = k^,T / (6nna) (/x being the fluid viscosity), we define 
two rescaled coupling diffusivities, 



A ± (C) 



D± {Q - Ds 
(a/R)D 



£ = x/R, 



(2) 



which are experimentally measurable and, within our far- 
field approximation, parameter-free (apart from the un- 
known geometrical factor 7). 

We now substitute the known solution for the Oseen 
tensor at the center of a cylindrical tube jl4|, |l5| to obtain 

00 

A ± (0 = ±(3/4) J2K cos(/U) + b n sin(/3 n 0]e- Q " ? 



±3/(2?), 



(3) 



" [ ±(3/4)[a 1 cos(/3 1 + &isin(/3 1 0]e- Ql «, £»1 

In Eq. (| + i[3 n are the complex roots 

of the equation u[Jq(u) + Ji(u)] = 2J (u)Ji(u), and 
(a n + ib n ) = 2{ir[2J 1 (u n )Y (u n ) - u n (J (u n )Y (u n ) + 
Ji( u n)Yi(u n ))]-u n }/J^(u n ) [J k and Y k being the Bessel 
functions of the first and second (Neumann) kind] . The 
coupling as a function of distance is depicted in Fig. [|. 
For ^ < 1 the particles are insensitive to the walls and 
the coupling approaches the algebraic ~ l/£ dependence 
as in an unbounded fluid frll| |. Yet, for £ > 1 the confined 
geometry becomes manifest; the sum is dominated by its 
first term and the interaction decays exponentially with 
distance. (The coefficients of this term are a\ ~ —0.0370, 
61 ~ 13.8, ai ~ 4.47, and ft ~ 1.47.) 

The sharp screening beyond the confinement length 
is unique to the linear geometry and is caused by the 
boundary conditions at the channel walls. We may visu- 
alize the boundary conditions as replaced by an infinite 




FIG. 3: Rescaled coupling diffusivity A + (= — A~) as a func- 
tion of rescaled inter-particle distance ^ — x/R [Eq. solid 
curve]. At very small distances the curve approaches the al- 
gebraic ~ l/£ dependence as in an unbounded fluid (dotted 
line) . At large distances it decays exponentially (dashed line) . 
For £ > 2.14 the coupling becomes negative (inset). 



series of image forces transverse to the channel axis, ac- 
companying the moving particle fl2]| . At distances larger 
than the channel width the images cancel the effect of the 
actual force. Small oscillations on top of the exponential 
decay make the coupling change sign at £ ~ 2.14 (Fig. ||, 
inset). This is a very peculiar effect, implying that for 
large distances the particles inhibit each other's motion 
rather than aid it. Physically, this arises from flow rolls 
that form along the channel pq| . Unfortunately, as seen 
in Fig. ^, even if the effect still exists in our partially open 
channel, it is far too small to be currently detectable jl9| . 

In Fig. ^ we have replotted the experimental data of 
Fig. [| scaled according to Eq. (||), with 7 = 2.44 for 
the square cross-section (samples 1 and 2) and 7 = 2.72 
for the rectangular one (sample 3). The data for dif- 
ferent channel widths, sphere sizes, and densities col- 
lapse onto two universal curves for the collective and 
relative diffusivities. The collapse confirms that the ob- 
served coupling is well described by a two-body hydro- 
dynamic interaction within a first-order approximation 
in a/w. Moreover, considering our crude approxima- 
tions, the agreement between the universal curve and 
Eq. (||) is remarkable. The stokeslet approximation thus 
provides a surprisingly reasonable description of the hy- 
drodynamic interaction even for relatively large parti- 
cles and short inter-particle distances. We can further 
use the fitted 7 to calculate the expected sclf-diffusivity 
which, at the same level of approximation pj ], is given 
by D s ~ A)(l-2.10444a/i?), R = 710/2. This yields for 
samples 1, 2 and 3D S ~ 0.16, 0.045 and 0.040 /im 2 /s, re- 
spectively. Although these values are not far off the mea- 
sured (low-density) ones — 0.11, 0.036 and 0.028 /im 2 /s — 
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the differences, compared to the good agreement in Fig. 
|], indicate that the self-diffusivity is more sensitive than 
the two-particle coupling to higher orders in a/w p0|. 
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FIG. 4: Rescaled coupling diffusivities A as a function of 
rescaled inter-particle distance £ = x/R for all the data of 
Fig. ^. The solid line is the theoretical curve for a particle in 
a cylindrical tube (Fig. ^|). The inset shows a cross-section of 
the sample cell with its theoretical effective cylinder. 



A small but systematic discrepancy between the cal- 
culation and experiment is seen in Fig. ^ as the distance 
becomes very small. This deviation naturally marks the 
breakdown of the x 3> a assumption. The observed trend 
can be understood by considering the extreme case of two 
particles in contact. The relative diffusivity then van- 
ishes, hence A - — » —D B / (Doa/R), whereas the collective 
one becomes equal to twice the self-diffusivity of a rigid 
"train" of two touching particles, 2D 2 - Since, evidently, 
D 2 < D a , we have A+ -► (2D 2 - D s )/(D a/R) < |A"|, 
in accord with the observed asymmetry. 

Another important consequence of the screened cou- 
pling is that concerted motion of neighboring particles 
becomes appreciable at a sharply defined density, com- 
parable to X/R) this dynamic clustering will be studied in 
more detail in a future publication. Our results suggest 
that a tube-stokeslet approximation [i.e., Eq. (0)] should 
give a good estimate of the hydrodynamic interaction in 
any linear geometry, provided that the particles are re- 
pelled from the walls and one has an estimate for the 
geometrical factor 7. Although we have discussed only 
freely diffusing particles, since mobilities are proportional 
to diffusivities, the same conclusions apply to driven mo- 
tion as well pH ], and should also bear upon motion and 
patterning in microfluidic systems p3] . 
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